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Abstract 

We studied nonequilibrium oscillations between left-handed electron neutrinos and 
nonthermalized sterile neutrinos in medium for resonant and nonresonant cases. The 
exact kinetic equations for neutrinos, written in terms of neutrino density matrix in 
momentum space were analyzed. The evolution of neutrino density matrix was numer- 
ically calculated. This allowed to study precisely the evolution of the neutrino number 
densities, energy spectrum distortion and the asymmetry between neutrinos and an- 
tineutrinos for each momentum mode. The influence of nonequilibrium oscillations 
on the primordial production of 4 He was calculated. Cosmological constraints on the 
oscillation parameters were obtained. 
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1 Nonequilibrium neutrino oscillations 



We discuss nonequilibrium oscillations between weak interacting electron neutrinos v e and 
sterile neutrinos v a for the case when v s do not thermalize till 2 MeV and oscillations become 
effective after u e decoupling. Oscillations of that type, but for the case of v a thermalizing 
before or around 2 MeV have been already discussed in literature [1-8]. We have provided a 
proper kinetic analysis of the neutrino evolution in terms of kinetic equations for the neutrino 
density matrix in momentum space. Such approach, though much more complicated, reveals 
some beautiful features of noneqilibrium oscillations, that cannot be catched otherwise. 
The assumptions of the model are the following: 

• Singlet neutrinos decouple much earlier than the active neutrinos do: Tjf > T^, as 
far as they do not participate in the ordinary weak interactions. 

Therefore, in later epochs their temperature and number densities are considerably less than 
those of the active neutrinos due to the subsequent annihilations and decays of particles 
that have additionally heated v e in equilibrium in comparison with the already decoupled 
v s - 

• We consider oscillations between the sterile v s and the active v e neutrinos, according 
to the Majorana&Dirac (MSzD) mixing scheme [9]. 

For simplicity in what follows we accept a simple M&D-mixing scheme with mixing present 
just in the electron sector z/j = Uu vi, I = e, s: f\ 

V\ = cv e + sv s 

Vi = SV e + Cf s , 

where v s denotes the sterile electron antineutrino, c = eos(i?) s = sm(i?) and i? is the 
mixing angle in the electron sector, the mass eigenstates v\ and V2 are Majorana particles 
with masses correspondingly m\ and m^. 

When transitions between active and sterile neutrinos are allowed, sterile neutrinos may 
not decouple from the plasma earlier enough in comparisson with the active ones, so that 
their densities may not differ enough for a noneqilibrium to be considerable, or they may 
regain their thermal equilibrium if already decoupled [11, 1-4]. The reactions of v e with 
the plasma are the source of thermalization for u s , because when oscillating into active 
neutrinos they may interact with the medium and thus thermalize. Therefore, 

• we assume that neutrino oscillations become effective after the decoupling of the active 
neutrinos, r osc > H for T < 2MeV. 

This puts constraint on the neutrino mass difference: 5m 2 < 1.3 x 10~ 7 eV 2 . 

• Finally we require that sterile neutrinos have not thermalized till 2 MeV when oscil- 
lations become effective. Q 

2 The transitions between different neutrino flavours were proved to have negligibie effect on the neutrino 
number densities and on nucieosynthesis because of the very siight deviation from equilibrium in that case 
Tf ~ T' f (/ is the flalour index) [10]. 

3 For the estimation of this rate we have used the pioneer work of Manohar [11] and also recent publications 
on the problem of thermalization of neutrinos by Barbieri&Dolgov [1], Kainulainen [2] and Enqvist et al. 
[4]- 
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This puts the following limit on the allowed range of oscillation parameters: sin{2'd)5m 2 < 

io- 7 e y 2 . 

As far as for this model the rates of expansion of the Universe, neutrino oscillations and 
neutrino interactions with the medium may be comparable, we have used kinetic equations 
for neutrinos accounting simultaneously for the participation of neutrinos into expansion, 
oscillations and interactions with the medium. Moreover, as far as nonequilibrium oscilla- 
tions are concerned, the density matrix of neutrinos considerably differs from its equilibrium 
form, and therefore one should work in terms of density matrix of neutrinos in momentum 
space [10, 12, 13]. ^ We have analyzed the evolution of nonequilibrium oscillating neutrinos 
by numerically integrating the kinetic equations for the density matrix in momentum space 
for the period after the decoupling of the electron neutrino till the freezing of neutron- 
proton ratio (n/p-ratio), i.e. for 2 MeV > T > 0.3 MeV. We considered both resonant 
5m? = m\ — m\ < and nonresonant 5m 2 > neutrino oscillations. 

Three main effects of neutrino nonequilibrium oscillations were revealed and precisely 
studied: 

(a) As far as oscillations become effective when the number densities of v e are much 
greater than those of u s , N Ue 3> N Vs , the oscillations tend to reestablish the statitstical 
eqilibrium between different oscillating species. As a result N Ve will decrease in comparison 
to its standard eqilibrium value. (This result of our study is in accordance with other 
publications concerning depletion of electron neutrinos [3-7]). 

(b) For the case of strongly noneqilibrium oscillations the distortion of the energy distri- 
bution of neutrinos may be considerable and must be accounted for properly^] This effect 
was not discussed in publications concerning active-sterile neutrino oscillations, and was 
thought to be negligible. In Ref. [12] it was first shown that for the case of v e <-> v s vacuum 
oscillations this effect is considerable and may be even greater than that of an additional 
neutrino species. In the present work we discuss this effect for neutrino oscillations in a 
medium. We have studied the evolution of the energy distortion using kinetic equations for 
neutrino density matrix in momentum space [10, 12, 13]. As far as these integro-differential 
equations are coupled nonlinear equations, analytic solution cannot be obtained without 
drastic assumptions. Therefore, we have numerically explored the problem. 

The evolution of the distortion is the following: First the low energy part of the spec- 
trum is distorted, and later on this distortion concerns neutrinos with higher and higher 
energies (Fig. 1). This behavior is natural, as far as neutrino oscillations affect first low 
energy neutrinos, r osc ~ 5m? /E v . The naive account of this effect by shifting the effective 
temperature and assuming the neutrino spectrum of equilibrium form gives wrong results 
for the case 5m? < 10~ 7 eV 2 . For bigger neutrino mass differences oscillations are fast 
enough and the naive account looks more acceptable, provided that u e have not decoupled. 

(c) Other interesting effect revealed by our approach is the generation of asymmetry 
between v e and their antiparticles. The possibility of an asymmetry generation due to 
CP-violating flavour oscillations was first proposed in Ref. [14]. Later estimations of an 
asymmetry due to CP-violating MSW resonant oscillations were provided [15]. Recently 
the problem of asymmetry was considered in connection with the exploration of the neutrino 

4 When neutrinos are in equilibrium their density matrix has its equilibrium form, namely = 
5ij exp(/i/T — E/T), so that one can work with particle densities instead of p. In an equilibrium back- 
ground, the introduction of oscillations slightly shifts p from its diagonal form, due to the extreme smallness 
of the neutrino mass in comparison with the characteristic temperatures and to the fact that equilibrium 
distribution of massless particles is not changed by the expansion [10]. 

This effect was discussed for a first time in Ref. [10], but unfortunately, as far as the case of flavour 
neutrino oscillations were considered and the energy distortion for that case was shown to be negligible, it 
was not paid the necessary attention it deserved. 
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propagation in the CP-odd plasma of the early Universe [1-5] and this type of asymmetry was 
obtained to be negligible. Recently it was shown in Ref. [16, 7], that asymmetry can grow to 
a considerable values for the case of great mass differences, 5m? > IQ^eV 2 was discussed. 
Our approach allows precise description of the asymmetry and its evolution, as far as 
working with the self consistent kinetic equations for neutrinos in momentum space enabled 
us to calculate the behavior of asymmetry at each momentum. So, the calculated result 
may differ considerably from the rough estimations made by working with neutrino mean 
energy and with the integrated quantities like particle densities and the energy densities. 

For the nonresonant case we obtained that the asymmetry is negligible. The asymmetry 
effect is noticeable only for the resonant case, as it was expected. Even when the asymmetry 
is assumed initially negligibly small (of the order of the baryon one), i.e. ~ 1CP 10 , it may 
be considerably amplified at resonant transition due to different interactions of neutrinos 
and antineutrinos with the CP-odd medium. The value of the asymmetry may increase 
by several orders of magnitude, oscillating and sign changing. This behavior may remain 
such long after the resonant transition. (This behavior needs very high accuracy and thus 
complicates strongly the numerical analysis.) Even in case when the value of the asymmetry 
does not become considerable enough to have some direct noticeable effect, on primordial 
nucleosynthesis for example, the asymmetry term at the resonant transition determines 
the evolution of the neutrino density matrix. (As far as the leading order terms, namely 
of oscillation effects and interaction with the medium, compensate each other at resonant 
transition.) It effectively suppresses the transitions of active neutrinos (antineutrinos) thus 
weakening the effect of neutrino depletion at resonance. Thus the evolution of neutrino 
may drastically differ from the case without asymmetry. For some model parameters this 
effect consists 20% of the previously discussed. We will discuss asymmetry evolution in 
more detail elsewhere. Here finally we want only to note, that linear approximation for the 
asymmetry in the resonant case may give very misleading results. 

As an illustration of the strength and importance of these effects we discuss their influ- 
ence on the primordial production of 4 He. Effect of oscillations on nucleosynthesis has been 
discussed in numerous publications [1-8, 12, 14]. However, a detail kinetic calculations of 
helium abundance for the case of nonequilibrium oscillations in medium has not been done. 
We describe it in the last section. Here we briefly summerize the results: 

The depletion of the electron neutrino number densities due to oscillations to sterile ones 
leads to an effective decrease in the weak processes rates, and thus to an increase of the 
freezing temperature of the n/p-ratio and corresponding overproduction of the primordially 
produced 4 He. The effect of the distortion of the energy distribution of neutrinos has two 
aspects. An average decrease of the energy of active neutrinos should lead to a decrease 
of the weak reactions rate , ~ E 2 and subsequently to an increase in the freezing 
temperature and the produced helium. On the other hand, there exists an energy thershold 
for the reaction v e + p — > n + e + . And in case when, due to oscillations, the energy 
of the relatively greater part of neutrinos becomes smaller than that threshold the n/p- 
freezing ratio decreases [17]. However, the numerical analysis showed that the latter effect 
is less noticeable compared with the previously described ones. The asymmetry calculations 
showed a slight predominance of neutrinos over antineutrinos, leading to decrease of helium. 
The greater effect of the asymmetry is, however, the relative increase of both neutrino and 
antineutrino particle densities which leads to a noticeable underproduction of helium (up 
to 10% relative decrease) Thus, the total result of nonequilibrium neutrino oscillations 
is overproduction of helium in comparison to the standard value. We have used the 4% 
relative increase in the primordially produced helium to obtain the exclusion region for the 
oscillation parameters. 
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2 The kinetics of noneqilibrium neutrino oscillations 



The kinetic equations for the density matrix of the noneqilibrium oscillating neutrinos in 
the primeval plasma of the Universe in the epoch previous to nucleosynthesis, i.e. consisting 
of photons, neutrinos, electrons, and small quantities of nucleons, have the form: 



dp{t) 
dt 



Hp 



dp® 
dp 



+ i [H , p(t)] + i [Hint, p(t)] + 0[H 



L int 



(1) 



where p is the momentum of electron neutrino and p is the density matrix of the massive 
Majorana neutrinos in momentum space. 

The first term in the equation describes the effect of expansion, the second is responsible 
for oscillations, the third accounts for forward neutrino scattering off the medium. H is 
the free neutrino hamiltonian: 



H Q 



P 



rrii 











p 2 + m| 



while Hint = ot V is the interaction hamiltonian, where ctij = U* e Uj e , V = Gf (±L — ^ / iv± w j , 
and in the interaction basis plays the role of induced squared mass for electron neutrinos: 



Q/M 2 



H 



LR 

int 




The first 'local' term in V accounts for charged- and neutral-current tree-level interactions 
with medium protons, neutrons, electrons and positrons, neutrinos and antineutrinos. It is 



proportional to the fermion asymmetry of the plasma L 
order of the baryon one 



L 



f 



N f -N 



f 



N B - N 



B 



Y^f Lf, which is assumed of the 



/3T 3 



The second 'nonlocal' term arises as an W/Z propagator effect, Q ~ E v T 4 . For the early 
Universe conditions both terms must be accounted for because although the second term is 
of the second power of Gf , the first term is proportional besides to the first power of Gf, to 
the very small value of the fermion asymmetry (of the order of the baryon asymmetry, i.e. 
1CP 10 in case B — L conservation is accepted) [18]. Moreover, the two terms have different 
temperature dependence and an interesting interplay between them during the cooling of 
the Universe is observed. 

The last term in the Eq. (|]) describes the weak interactions of neutrinos with the 
medium. For example, for the weak reactions of neutrinos with electrons and positrons 
e + e~ 



e'^v'i it has the form 



/ dQ(u, e + , e 
+ JdQ(e + ,v',e' + ) 



n,-n P + 



AA^ - \ [p, A^pA) 
> e _Bp'B^ -\{EtB, p} + r, 
'e+Cp'tf -|{ctC, p) n t 



where n stands for the number density of the interacting particles, A is the amplitude of 



the process e + e 



B is the amplitude of the process e Vj 



and C is the 



4 



amplitude of the process e + Vj — > e' + u[. They are expressed through the known amplitudes 

A e (e + e~ -» u e u e ), B e (e~v e ->• e~z/ e ) and C e (e+^ e -> e+z^ e ): 

A = a A e , B = a B e . C = a C e . 

We have analysed the evolution of the neutrino density matrix assumed that oscillations 
become noticeable after electron neutrinos decoupling. So, the neutrino kinetics down 
to 2 MeV does not differ from the standard case, i.e. electron neutrinos maintain their 
equilibrium distribution, while sterile neutrinos are absent. Then the last term in the 
kinetic equation can be neglected. The equation results into a set of coupled nonlinear 
integro-differential equations for the components of the density matrix. Analytical solution 
is not possible. [] 

We have numerically calculated the evolution of the neutrino density matrix for temper- 
ature interval [0.3,2.0] MeV, i.e. after neutrino decoupling and till the n/p freezing. The 
oscillation parameters range studied is 5m 2 € ±[10 -10 , 10 -7 ] eV 2 and § £ [0,7r/4]. The 
baryon asymmetry (3 was taken to be 3 x 10 -10 . 

As was already stated our numerical analysis showed that the nonequilibrium oscillations 
can considerably deplete the number densities of electron neutrinos (antineutrinos) and 
distort their energy spectrum. Neutrino-antineutrino asymmetry may grow at the resonant 
transition and may also effect considerably the evolution of neutrino ensembles. 



3 Nucleosynthesis with nonequilibrium oscillating neutrinos 

The kinetic equation describing the evolution of the neutron number density in momentum 
space n n for the case of oscillating neutrinos v e «-> v s reads: 

lt = H Pn%t +Jdn(e-,p,iy)\A(e-p^vn)\ 2 (n e -n p -n nPLL ) 
-J dU(e + ,p, 9)\A(e + n —> pv)\ 2 (n e +n n -n p p LL ). 

The first term describes the expansion and the second one - the processes e~+p <-> n-\-v e and 
p+i> e <-> e + +n, directly influencing the nucleon density. It differs from the standard scenario 
one only by the substitution of pll and pll instead of n e u q ~ exp(—E u /T). The neutrino 
and antineutrino density matrices may differ pn, ^ pll, contrary to the standard model, as 
a result of the different reactions with the CP-odd plasma of the prenucleosynthesis epoch. 
Number densities per unit volume are expressed as N = (27r) -3 / d 3 p n{p). 

We have numerically integrated the equation for the temperature range of interest T 6 
[0.3, 2.0] MeV for the range of oscillation parameters of our model. The results of the 
numerical integration are illustrated in Fig. 2, in comparison with the vacuum case and the 
standard nucleosynthesis without oscillations. It is obvious that the kinetic effects of the 
nonequilibrium oscillations on nucleosynthesis are considerable. The next figure presents 
the dependence of the primordially produced helium on mixing angle (Fig. 3). 

From numerical integration for different oscillation parameters we have obtained con- 
stant 4 He contours. The Fig. 4 gives the cosmologically excluded regions of oscillation 
parameters. It corresponds to primordial abundance of helium Y p = 0.245, which gives 4% 
overproduction of helium in comparison with the observational value. 

For the cases when the energy distortion and asymmetry are considerable we have ob- 
tained an order of magnitude stronger constraints than the cited in literature [4-7]. There- 
fore, in conclusion we would like to stress once again, that in case of nonequilibrium neutrino 

6 For the case of vacuum neutrino oscillations this equation was analytically solved and the evolution of 
density matrix was given explicitly in Ref. [12]. 
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oscillations working with the exact kinetic equations for the density matrix of neutrinos in 
momentum space is necessary. 
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Figure Captions 



Figure 1: The figure shows the energy distortion of active neutrinos x 2 pll(x), where 
x = Ey/T, for the case of nonequilibrium neutrino oscillations, 5m? = — 10 -8 , $ = ir/8 at 
different temperatures: T = 1 MeV (a), T = 0.7 MeV (b), and T = 0.5 MeV (c). 

Figure 2: The curves represent the evolution of the neutron number density relative to 
nucleons X n {t) = N n (t)/(N P + N n ) for the nucleosynthesis model with vacuum nonequilib- 
rium oscillations and for the case of nonequilibrium oscillations in medium, 5m? = — 10~ 8 , 
•d = 7r/8. For comparison the curve corresponding to the standard nucleosynthesis model is 
shown. 



Figure 3: The figure illustrates the dependence of the neutron number density relative 
to nucleons X n = N n /(N p + N n ) on the mixing angle for 5m? = — 10 -8 . 

Figure 4: Exclusion regions for oscillation parameters are shown for the case of resonant 
5m? < and nonresonant 5m? > neutrino oscillations. The curves correspond to helium 
abundance Y p = 0.245. 



